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Abstract 



We study J /if) suppression in AB collisions assuming that the charmonium 
states evolve from small, color transparent configurations. Their interac- 
tion with nucleons and nonequilibrated, secondary hadrons is simulated us- 
ing the microscopic model UrQMD. The Drell-Yan lepton pair yield and the 
J/^/Drell-Yan ratio are calculated as a function of the neutral transverse en- 
ergy in Pb+Pb collisions at 160 GeV and found to be in reasonable agreement 
with existing data. 
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I. INTRODUCTION 



Experimental data on the AB and Ex dependence of J/^-meson production from the 
CERN SPS (l|-|4| exhibit tantalizing evidence for J/ip suppression when compared to hard 
QCD production. Do the J ftp data in Pb+Pb-collisions indicate the creation of a deconfined 
phase of strongly interacting matter, i.e. a quark-gluon plasma (QGP) 0,0? Such QGP 
scenarios generally rely on the experimental observation of deviations from the predictions 
of models of hadronic suppression. Simple analytical models of J/ip suppression have used 
nuclear absorption alone || or absorption by nucleons and comoving secondaries j7| to fit 
all but the Pb+Pb data. These models generally assume a single, fixed nucleon absorption 
cross section for all charmonium states and that the comover interaction rate, assuming a 
Bjorken scaling expansion, is essentially thermal. (For theoretical reviews, see Refs. []8| [Toll . ^)F1 

Charmonium final state interactions with mesons and baryons are simulated using the 
Ultrarelativistic Quantum Molecular Dynamics, UrQMD ||15|| , a microscopic hadronic trans- 
port model. The charmonium nucleon cross sections employed here are calculated in a 
nonrelativistic potential model [16|] since at SPS energies, charmonium-nucleon interactions 
are predominantly nonperturbative. The charmonium precursors are not eigenstates of the 
QCD Hamiltonian, therefore their effective sizes and interaction cross sections may vary from 
their production to the final state formation flTEfl . Thus the charmonium absorption cross 
sections are assumed to expand linearly with time until their asymptotic value is reached. 
Since the produced particles in the collision have a non-thermal distribution [I7|, comover 
interactions can rather effectively dissociate the J/ip. Because of abundant high mass res- 
onances, most of the meson induced J/ip dissociation processes are exothermic ||17|| . They 



account for an important fraction of the observed suppression, larger in Pb+Pb collisions 
than the nuclear absorption alone. 



II. THE MODEL 



We apply perturbative QCD to the production of charmonium states by simulating 
nucleus-nucleus collisions in the impulse approximation. The nuclear dependence of parton 
distribution functions is neglected. The resulting space-time distribution of charmonium 
production points is inserted into the evolving hadronic environment calculated with UrQMD 
fP7f since the rare quarkonium production processes are small perturbations on the heavy 
ion collision^. Our model is thus designed to account for partonic and hadronic aspects of 
the charmonium dynamics. The space-time distributions of hard and soft processes, i.e. 



1 Recently, there have been several attempts to build models of charmonium production and ab- 



sorption by means of microscopic hadronic transport simulations [11-14]. Based on rather different 
model treatments of the charmonium dynamics and the interconnection of hard and soft processes 
it is claimed in all these studies that conventional hadronic scenarios are consistent with the Pb+Pb 
data. 

2 Since we consider only exclusive charmonium production, the QCD factorization theorem is 
inapplicable. 
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in particular charmonium production and absorption, may overlap. Thus, in our model, a 
charmonium state produced in a hard process can be dissociated by the interaction with a 
comoving hadron before this state leaves the nuclear environment and before all other hard 
production processes are completed (in contrast to Ref. JL4[]). The probability for such an 
event, however, is reduced due to the initially small absorption cross sections and the finite 
formation times of comoving mesons according to the string model (see below). To avoid 
double counting, interactions of cc states with produced hadrons in individual NN collisions 
are excluded. In AB collisions, all hard NN collisions can contribute to the cc production 
while all soft iViV collisions can contribute to the hadronic environment in which the cc state 
may be dissociated. The error imposed by this concept — inherent to all microscopic and 
analytical models of comover absorption — is estimated to be very small. 

The charmonium states are distributed according to their assumed production probability 
times their decay probability to J/i/j's. Thus 40% of the final states are x's, 55% are J/ip's, 
and 5% are ^'s [|TH|. According to the spin degeneracy, 1/3 of the x' s are Xcio states and 



2/3 are Xcii states. Their momenta are assigned according to the parametrization |T9 



The rescattering cross sections for X(cc) + B, assuming B = N, are taken from Ref. [16 
a(J/ipN) = 3.6 mb, a{ijj'N) = 20 mb, a(XcioN) = 6.8 mb, and a(xdiN) = 15.9 mb. 
Charmonium- meson cross sections (X(cc) + 7T, X(cc) + p, etc.) are reduced by a factor 
of 2/3 from the corresponding baryon values. All baryon and meson collisions above the 
respective dissociation thresholds are assumed to break up the charmonium state. Universal 
and energy independent cross sections are employed, ignoring any charmonium-meson reso- 
nances, perhaps too crude an assumption. From phase space arguments one can infer that 
the J/ip dissociation cross sections with 7r's should be suppressed close to threshold while 
it should be enhanced for exothermic channels. However, we have found that, during the 
initial stage of the J/ip comover collisions, the average interaction energy, < E >k, 5 GeV, 
is far above threshold |17| . There are no calculations of J / ip dissociation cross sections with 



mesons other than 7r's and p's JTI|. In a thermal comover scenario the density of heavier 
mesons is suppressed by the Boltzmann factor. According to the UrQMD simulations, how- 
ever, these mesons dominate the comover absorption |I7j]. In a forthcoming paper [^0J the 



influence of the energy dependence of the comover interactions will be studied further. 

The cross sections correspond to the geometrical transverse radii r % T = of the char- 
monium states. We use a % to estimate the respective formation times r F of the charmonium 
states by choosing t 1 f = r l T /c. During these formation times the cross sections increase 



linearly with t ||16| , starting from at t = 0. 

Here it is important that we also take into account the formation time of comoving 
mesons (on average, tf ~ 1 fm/c). Particles produced by string fragmentation are not 
allowed to interact with other hadrons - in particular with a charmonium state - within 
their formation time. However, leading hadrons are allowed to interact with a reduced cross 
section even within their formation time. The reduction factor is 1/2 for mesons which 
contain a leading constituent quark from an incident nucleon and 2/3 for baryons which 
contain a leading diquark. 
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For this study, we have slightly modified the angular distributions of meson-baryon inter- 
action in the UrQMD 1.0 model since their strong forward peak underpredicts the total pro- 
duced transverse energy The model now reproduces the Et spectra in S(200 GeV)+Au 
and Pb(160 GeV)+Pb collisions measured by NA35 and NA49, respectively Neither 
the amount of baryon stopping nor the rapidity distribution of negatively charged parti- 
cles which have been shown to agree with experimental pp and AB interactions 
significantly affected by this modification. 
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are 



III. RESULTS AND DISCUSSION 

Figure [l] shows the calculated number of Drell-Yan muon pairs, proportional to the num- 
ber of hard collisions, in Pb+Pb collisions as a function of the produced neutral transverse 



energy within 1.1 < r\ < 2.3. The NA50 data |23] have been included with the abscissa 
rescaled to reflect the latest change in data M which indicates an w 20 % shift in the abso- 
lute Et scale from previous publications [j3|,p3"| . We are aware that the new analysis does not 
imply a simple overall rescaling of the old data points. However, in order to reasonably com- 
pare the gross features of the experimental dimuon Et spectrum with our model calculation, 
we have multiplied all -EV-values of the data by 0.8. This factor was obtained by comparing 
the E T — E Z dc contour from Quark Matter '97 [^3| and the Moriond '98 [|J proceedings. 
The agreement between the model and the rescaled data is to be expected since the NA49 
Et distribution |MJ is described correctly |22[0. The additional simulation is a simple and 



well understood model of hard scattering processes in nucleus-nucleus collisions. 

Figure |2] shows the J/ ip production cross section according to UrQMD calculations for 
several projectile-target combinations (p(450 GeV)+C, p(200 GeV)+Cu, p(200 GeV)+W, 
p(200 GeV)+U, S(200 GeV)+U, and Pb(160 GeV)+Pb) in comparison to experimental 
data [[3]. The results of the calculations are normalized to the experimental cross section 
in p(200 GeV)+p interactions. The 450 GeV and 160 GeV simulations are rescaled to 
Pi a b = 200 GeV with the parametrization of Ref. [PSfl , as done by NA50 Considering 
only nuclear dissociation results in a far smaller J/ ip suppression than seen in the data, not 
only for Pb+Pb collisions but also for S+U and even pA reactions. Note that the system- 
atics of nuclear absorption shown in Fig. |2] does not reflect a universal straight line as in 
Glauber calculations with constant absorption cross sections. By taking the nonequilibrium 
charmonium- meson interactions into account, good agreement with the data is obtained. 
However, a strong dependence on parameters such as the charmonium and comover forma- 
tion times and the dissociation cross sections remains to be studied in detail j2(J . 

Due to the linear expansion of the charmonium cross sections with time, the J/ip and 



3 However, the agreement between the model and the NA38 experiment becomes poor for S+U 
collisions. The UrQMD calculation appears to overestimate the neutral transverse energy by 
about 25% in the range 1.7 < r) < 4.1 flj ^l although the calculated Et spectrum of the similar 
S+Au system agrees well with the NA35 data [24]. The UrQMD calculation thus indicates an 
inconsistency between the Et measurements by NA35, NA49 and NA50 on the one hand and 
NA38 on the other |2l. 
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ip' cross sections are similar in the very early stages, leading to a weak A- dependence of 
the ip' / J/ip ratio in pA collisions at central rapidities, apparently consistent with the data. 
However, for a deeper understanding of this ratio quantum interference effects p7[ as well 
as refeeding processes, n J / ip — > ip'n [p8|p9| , must be considered. The if)' '/ J / ip ratio will be 
studied further later pOf . 

Figure ^| shows the J / ip to Drell-Yan ratio as a function of Et for Pb+Pb interactions 
at 160 GeV compared to the NA50 data |§. The normalization of B^a( J/ip)/cr(DY) = 46 
in pp interactions at 200 GeV has been fit to S+U data within a geometrical model ||. The 
application of this value to our analysis is not arbitrary: the model of Ref. |5| renders the 
identical ^-integrated J/ip survival probability, S = 0.49, as the UrQMD calculation for 
this system. An additional factor of 1.25 [f§ has been applied to the Pb+Pb calculation in 
order to account for the lower energy, 160 GeV, since the J/ip and Drell-Yan cross sections 
have different energy and isospin dependencies. The gross features of the E T dependence 
of the J/ip to Drell-Yan ratio are reasonably well described by the model calculation. No 
discontinuities in the shape of the ratio as a function of Et are predicted by the simulation. 



IV. CONCLUSION 

We have examined charmonium production and absorption processes in pA and AB 
collisions at SPS energies. The microscopic simulation of hard processes in the impulse 
approximation and the hadronic transport description of AB collisions with the UrQMD 
model simultaneously provide reasonable Et dependencies of the Drell-Yan rates as well as 
baryon and meson rapidity distributions. We have modelled J/ip absorption according to 



the scenario described in Ref. | 16[] . Cross sections evolving from color transparent small 
configurations to asymptotic states derived from quantum diffusion and the dynamical x 
polarization (color filtering) are taken into account. 

The calculated J/ip production cross sections for minimum bias pA, S+U and Pb+Pb 
collisions agree with experiment. Dissociation by nonequilibrium comovers accounts for 
about half of the total absorption in S+U and Pb+Pb reactions. The contribution of the 
interaction with comovers in pA reactions is small but not negligible. The suppression 
of charmonium states is sensitive to the comover momentum distributions. The effective 
dissociation by comovers seems to indicate a nonequilibrated hadronic environment. The 
observed Et dependence of the J/ip to Drell-Yan ratio in Pb+Pb collisions is reproduced 
by the model. The calculated result is smooth, without abrupt discontinuities, in agree- 



ment with new high statistics data |3(J. We conclude that within our model, the data on 



charmonium cross sections at the SPS can be explained without invoking exotic mechanisms. 
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FIG. 1. Number of Drell-Yan pairs in Pb+Pb interactions as a function of the neutral trans- 
verse energy within 1.1 < rj < 2.3. The calculation is normalized to the data. Shown is the 
UrQMD result and experimental data from NA50 [^3| with the Er of the data rescaled by 0.8. 
The modification is motivated by a comparison of the recently published Et — Ez dc contour plot 
|Q] with the previously published analysis ]3|,p3|]. 
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FIG. 2. J/^-production cross sections times dimuon branching ratio in the kinematical domain 
< Hem < 1 and | cos 9cs I < 0.5, and rescaled, if necessary, to pi a b = 200 GeV as a function of 
AB. The data (open triangles) are from ||. Open circles denote the production cross sections if 
only nuclear absorption is considered. 
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FIG. 3. The ratio of J /ip to Drell-Yan production as a function of Et for Pb+Pb at 160 GeV. 
The experimental data are from Ref. The normalization factor, from, pp interactions at 200 GeV, 
B )j,tj,o (J /ip) I 'a (DY) = 46 is taken from Ref. ]5). This value, however, has been indirectly determined 
in the framework of a different model. An additional factor of 1.25 || has been applied to the 
Pb+Pb calculation in order to account for the lower energy. Note that no scaling factor has been 
applied to the x-axis for either the calculations or the data. 
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